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A chemical potential difference between the legs of a two-leg ladder is found to be harmful for 
Cooper pairing. The instability of superconductivity in such systems is analyzed by compairing 
results of various analytical and numerical methods. Within a strong coupling approach for the t-J 
model, supplemented by exact numerical diagonalization, hole binding is found unstable beyond a 
finite, critical chemical potential difference. The spinon-holon mean field theory for the t-J model 
shows a clear reduction of the the BCS gaps upon increasing the chemical potential difference leading 
to a breakdown of superconductivity. Based on a renormalization group approach and Abelian 
bosonization, the doping dependent phase diagram for the weakly interacting Hubbard model with 
different chemical potentials was determined. 
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I. INTRODUCTION 

Doped spin liquids have been an important sub- 
ject of condensed matter research for the last two 
decades, mainly due to their possible relevance 
to the (cuprate) high-temperature superconductors 
(HTSCs)i^ Although it is still not understood how high- 
temperature superconductivity emerges from an antifer- 
romagnetic Mott-insulator upon carrier doping, there is 
broad consensus that an intermediate pseudogap phase 
plays a crucial role in understanding both the exotic nor- 
mal and superconducting properties of these materials. 3 

Many proposals concerning the nature of the pseudo- 
gap phase have emerged. One candidate for the pseudo- 
gap phase is the resonating valence bond state (RVB), 
which describes strongly fluctuating short-ranged spin 
singlets. 1 While the relevance of this state to the quasi- 
two-dimensional HTSCs is still under debate, it is re- 
alized in specially designed lattice structures^ Of par- 
ticular interest are systems with ladder-shaped crystal 
structures, which are realized in various transition metal 
oxide compoundsAii 

Undoped, these ladder-systems are Mott-insulators 
and well described by quantum spin-1/2 Heisenberg mod- 
els. For the two-leg ladder the ground state is an RVB- 
like state displaying short-ranged spin-singlet correla- 
tions with spin-singlet dimers on the rungs dominating 
over those along the legs. This constitutes a spin liquid 
with a finite energy gap to the lowest spin excitations. 
Furthermore, it was proposed that such a system would 
exhibit superconductivity upon hole-doping£& Various 
theoretical approaches confirmed a strong tendency to- 
wards formation of Cooper pairs with phase proper- 
ties reminding of the d :E 2_„2-wave channel of the two- 
dimensional HTSCsi 9 ! 10 ! 11 ! 1 ^ 

The theoretical proposals were followed by an intense 
material research attempting to dope holes into a vari- 
ety of known insulating copper-oxide compounds display- 
ing ladder structures^ Materials under consideration are 
SrCu203, CaCu203, LaCu02.5, and SruCu^O/ti. Dop- 
ing these compounds is intrinsically difficult because of 
chemical instabilities, and carrier localization effects may 
inhibit the desired metallic behavior. Nevertheless, the 



search for superconductivity has been successful in the 
compound Sri4_ x Ca x Cu.2404i which contains layers of 
ladders alternating with layers of single chains^ In this 
material T c rises up to about 12 K under high pressure, 
which apparently leads to a transfer of charge carriers 
from the chains onto the ladders. However, a detailed 
understanding of this system under high pressure has not 
been reached yet. A well-controlled and less invasive way 
of doping a ladder compound is thus highly desirable. 

Hole and electron doping by means of field effect de- 
vices induces mobile charge carriers into the originally in- 
sulating material using a large gate voltage. This method 
would be ideal for doping quasi one- and two-dimensional 
systems, since the induced charge is confined to the outer- 
most layer of the compound, closest to the gate. 

An alternative technique of tunable doping has been 
achieved using heterostructures of layered materials 
such high-T c -cuprates combined with ferroelectrics like 
Pb(Zr,Ti)Oj 15 i 16 i 17 . These non-volatile techniques of 
tunable doping may allow for a detailed comparison be- 
tween experiment and theory in these low-dimensional 
structures. 

The most natural choice of a ladder system for this 
type of doping controll is a film in which the ladder planes 
lie parallel to the gate or ferroelectric, so that the carriers 
enter the ladders uniformly. However, in the compound 
LaCu02.s the ladders are not parallel to each other, but 
exhibit a staggering^ Consequently, in a field effect de- 
vice the ladders would be inhomogeneously doped in the 
sense that the chemical potential on the two legs would be 
different. Similarly, a variable orientation of the dipolar 
moments of the ferroelectric can lead to inhomogeneous 
doping. 

In this paper we analyze the evolution of the supercon- 
ducting state under such doping circumstances. Our an- 
alytical and numerical analysis shows that non-uniform 
doping is harmful to the superconducting state of the 
two-leg ladder. While for small A/i the ladder remains 
superconducting, the pairing is suppressed upon increas- 
ing A/i, and depending on the doping level new phases 
with reduced, and eventually without superconductivity 
appear. 

The paper is organized according to the used ana- 
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lytical and numerical methods. In this way we study 
various aspects of the problem within different approx- 
imative schemes. In the following section a qualitative 
argument for the limited stability of the superconduct- 
ing state upon inhomogeneous doping is presented. In 
Sec. Ill, a discussion of numerical exact diagonalization 
results is given for the charge correlations of two holes 
in ladders with up to 22 sites. Then, in Sec. IV we con- 
sider a mean field treatment of the t-J model based on 
the spinon-holon decoupling scheme. In Sec. V, we apply 
renormalization group (RG) and bosonization methods 
to derive the phase diagram of the weakly interacting 
Hubbard model in the inhomogeneously doped case. We 
conclude in Sec. VI and draw a unified picture of the be- 
havior of inhomogeneously doped two-leg ladder systems 
by combining the results from the earlier sections. 



II. STRONG RUNG COUPLING LIMIT 

The influence of a difference in the chemical poten- 
tial on the pairing state of the two-leg ladder can be 
illustrated by a simple qualitative argument for the t-J 
model. Consider the two-leg ladder with electrons mov- 
ing along the legs and rungs with hopping matrix ele- 
ments t and t' respectively, and nearest-neighbor spin 
exchange with exchange constants J and J' . The Hamil- 
tonian of the t-J model then reads 

H = -tY,V(c} as c 3+ha , s +h.c.)T 

j,a,s 
3,0 

■ • / 'X! ::S '- 1 ' S ' J 7".- 

3 

~ ^ Ma n 3,a,s- (1) 

3,o,s 

The operator c t a s (cj. a . s ) creates (annihilates) an elec- 
tron with spin s on site (j, a), where j labels the rungs, 
and a = 1, 2 the legs. The electron number operators are 
defined as n J} a }S = Cj as Cj, a , s , and n^ a = Yl s n j,a,s- The 
spin operators are 

^j.a = 2 ^ ^ C \.a.s a ss' c j,a,s' (2) 

where er a , a — 1,2,3 are Pauli matrices. The constraint 
of excluded double occupancy is enforced by the projec- 
tion operator 

V = - n 3,o,f n 3^,l)- ( 3 ) 

3,0, 



In the last line of Eq. QJ different chemical potentials 
on the two legs, (x a , a = 1, 2 describe an inhomogeneous 
doping of the system. Throughout this paper we assume 
A/i = /xi — i-J.2 > 0, and refer to the leg with a = 1 
(a = 2) as the upper (lower) leg. Furthermore, the overall 
doping concentration 6 = 1 — n fixes the average chemical 
potential p, — (/ii + M2)/2- 

The phases of the t-J model on the two-leg ladder for 
A/i = are well characterized^ At half-filling (5 = 0) 
with one electron per site the ladder is a (Mott-) insu- 
lating spin liquid. Upon removal of electrons, i.e. doping 
of holes, mobile carriers appear, resulting in a Luther- 
Emery liquid with gapless charge modes and gapped 
spin excitations. Furthermore, the gapless charge mode 
exhibits dominant superconducting correlations with d- 
wave-like phase structure. 

We now discuss the effect of inhomogeneous doping on 
this superconducting state, described by A/x > 0. For 
many aspects of ladder systems the limit of strong rung- 
coupling gives useful insights into their basic properties. 
Therefore we first consider the Hamiltonian (1) in the 
limit J', t' 3> J, t. Neglecting the coupling along the legs 
entirely, the undoped system corresponds to a chain of 
decoupled rungs, and the groundstate becomes a product 
state of dimer spin-singlets on the rungs. Furthermore, 
the lowest spin excitation corresponds to exciting one 
rung-singlet to a triplet, at an energy expense of J'. The 
superconducting state, i.e. Cooper pairing, in the doped 
spin liquid is inferred from the fact that two doped holes 
rather reside on a single rung rather than to separate 
onto two rungs. This is the case if the dominant energy 
scale is the spin exchange interaction, J' . Then the cost 
of breaking two spin singlets is larger than the gain of 
kinetic energy from separating the two holes. Namely, 
for two holes on a single rung the energy is 

E 2h = J' - 2/2, (4) 

while for a single hole 

E lh = J'-fl-^ ^A(t'f + Afj, 2 . (5) 
Pairing on a rung is favored, if 

2E lh - E 2h = J' - v/4(f) 2 + A M 2 > 0, (6) 

which in the uniformly doped case (A/i = 0) leads to the 
condition J' > 2i' for pairing. Obviously, a finite value 
of A/i weakens the pairing by reducing the above energy 
gain. 

This simple observation of depairing under non- 
uniform doping is confirmed by more sophisticated ap- 
proaches, as those considered in the following sections. 

III. EXACT DIAGONALIZATION 

To extend the discussion of the two-hole problem 
beyond the strong coupling limit we performed ex- 
act diagonalizations of finite systems, using the Lanc- 
zos algorithmi^S We considered the Hamiltonian Q at 
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isotropic coupling (t — t' , J — J'), and studied the half- 
filled system doped with two holes, using periodic bound- 
ary conditions. We studied systems of different length, 
L, containing 8 to f 1 rungs, and furthermore considered 
different values of J/t. 

Consistent with the strong coupling argument of Sec. 
II, the hole bound state is found to be unstable beyond a 
critical value of A/i > A/i c . Furthermore for the range of 
parameters considered here (0.4 < J/t < 0.8), this crit- 
ical value is A/i c rj J' . This indicates that the physics 
of the system is quite well captured by the strong rung- 
coupling limit, with J' being the dominant energy scale 
for pairing. 
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FIG. 1: Rung hole-hole correlation-function, (n (0)n (j)) , in 
the ground state of the t-J model for an inhomogeneously 
doped two-leg ladder at selected values of the chemical po- 
tential difference, A/i. The finite ladder has L — 10 rungs, 
and furthermore t — t', and J = J' = 0.5t. 



The behavior of the holes under non-uniform doping 
can be analysed using the hole-hole correlation func- 
tion. Denoting the hole number operator on rung j 
by n'(j) — 2 — rij } i — %\2, the rung hole-hole corre- 
lation function is defined as the ground state expec- 
tation value (n'(O)n'(j)), for the rung-rung separation 
j = 0,1, ■ This correlation function is shown 
in Fig. \I\for a ladder with 10 rungs for J/t = 0.5, and 
at selected values of the chemical potential difference, 
A/i. The behavior of the correlation function changes 
abruptly between A/z/i = 0.5 and A/i/t = 0.6. For small 
values of Afi/t < 0.5, we find maximal rung hole-hole cor- 
relations between neighboring rungs, and a strong decay 
of the correlation function at larger distances. For val- 
ues of A/i/t > 0.6, the value of the correlation function 
is increasing with distance, and has a maximum value 
at the maximal possible rung-rung separation. Further- 
more, it almost vanishes on the same rung. This clearly 
indicates the destruction of the hole-hole bound state for 
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FIG. 2: Characteristic hole-hole separation length, £, in the 
ground state of the t-J model for an inhomogeneously doped 
two-leg ladder as a function of the difference in the chemical 
potential, A/i, for L — 8 (circles), and L — 10 (squares). 
Furthermore, t' = t, and J — J' — 0.5t. 



A/i/t > 0.6, where the system consists of two holes on 
the lower leg, favored by a lower chemical potential. 
The sudden change in the behavior of the correlation 
function is due to a level crossing at A/i c in this system. 
While in the bound regime the ground state has zero to- 
tal momentum along the legs, the lowest energy state in 
the unbound regime has a finite value of the total mo- 
mentum. The correlation functions for A/i/i = 0.2 and 
A/i/t = 0.05 are almost identical, reflecting the robust- 
ness of the bound state against small doping inhomogeni- 
ties. Indeed, it can be shown that the doping asymmetry 
A/x does not have any effect in first order perturbation 
theory. In the regime of unbound holes the correlation 
function is again insensitive to changes in A/i, since for 
large A/i the holes are almost exclusively located on the 
lower leg, and therefore increasing A/i merely results in 
an overall energy shift. 

The drastic change in the hole-hole correlation function 
as a result of the transition between the bound and un- 
bound regime is also reflected in the characteristic hole- 
hole separation length, £, defined by 



E • 



(7) 



where N is a normalization factor given by N — 
^2 An' (Q)n! (j)) . In Fig. [3 this separation length is plot- 
ted as a function of A/i for two different systems with an 
even number of rungs, L = 8 and L = 10^i For both sys- 
tems the characteristic hole-hole separation length jumps 
discontinuously at a critical value of A/i c w J' . Further- 
more, in the bound regime, i.efor A/i < A/i c , finite size 
effects in this quantity are rather small already for the 
system sizes considered here. In this regime the holes are 
bound, and the rung hole-hole correlation decays expo- 
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nentially. Therefore, the bound pair wave function ex- 
tends over only a few rungs. In the unbound regime 
however, the two holes tend to separate onto the lower 
leg, and therefore £ grows with increasing system size. 

The above numerical analysis demonstrates that non- 
uniform doping, by confining the mobile carriers onto one 
of the legs, is harmful to pairing. Furthermore, the inter- 
leg exchange interaction plays the important role of stabi- 
lizing the bound hole pair state. While the finite ladders 
considered in this section may be viewed as systems with 
a doping concentration of roughly 8 ~ 0.1, different ap- 
proaches are needed to analyze the finite-doping regime 
beyond the two-hole problem. These will be presented in 
the following sections. 



IV. MEAN FIELD ANALYSIS 

In this section we extend our analysis of the t- J-model 
by considering a mean field approximation based on the 
spinon-holon-decoupling scheme. We follow a similar ap- 
proach as various previous studies on ladders as well as 
two-dimensional systemsii2i22i22i2a 



A. Spinon-Holon Decomposition 

The non holonomic local constraint cj a s Cj, a , s < 
1 is one of the main difficulties in treating the t-J- 
model. The slave-boson formalism provides a possibil- 
ity to take this constraint into account. Introducing 
fermionic spinon operators / and bosonic holon opera- 
tors b, the electron creation and annihilation operators 
can be expressed as 

c la,s = /j,«,Aa> and 'J— = &],«/i,a,«, (8) 
leading to the holonomic constraint 

£/!«,./*»,. +&i a. = (9) 



where the Lagrange multipliers A ja , a = 1,2 have been 
introduced to enforce the local constraint JJjJ. In the in- 
teraction part, the density-density terms n^a^j',^ are 
omitted. Within the following mean field treatment, this 
term would destroy the local SU (2) gauge symmetry of 
the spinon representation at half-filling. 2 - This symme- 
try corresponds to a local unitary rotation of the spinor 
(/i.Siti/ial) leaving the spinon spectrum invariant^ 
This symmetry appears naturally in the large-JJ Hub- 
bard model 2 - and is considered to be essential for various 
aspects of the weakly doped t- J-modeli^ 2 * 2 ^ Therefore, 
we will keep only the spin exchange part of the inter- 
action which conserves this symmetry in the mean field 
approximation^ The last term in Eq. (|ll)fl takes the dif- 
ferent chemical potentials on the two legs into account. 

To proceed we decouple the terms which are not single 
particle terms in the Hamiltonian (|10fl by introducing the 
following mean fields^ 



Xj-a;j',a' — „ ^^Sfj.a.s fj',a',s), 
s 

Bj,a;j',a' — (bj,a,bj, a t) , 
Aj,a;f,a' — (fj,a,lfj',a',y- 



(11) 



In the following the mean fields along the legs are labeled 
with the indices 1 and 2, and the mean field on the rung 
with the index 3, e.g. for \ : 

Xa = ^ /i+W> a = 1,2 (12) 

S 

X3 = o y^(/).l.s /j.2,a)- 
s 

This convention of labeling the bond (j, a;j', a') also ap- 
plies to the mean fields B, and A. Finally the doping 
level is fixed by the condition 



The Hamiltonian can be expressed in terms of this 
new operators as 

j,s 

+ J ^2 S j,a ■ S j+l,a 

j 

■ £ l> :>>],, ''j.- (10) 



(13) 



The Lagrange multipliers A JQ are kept uniform on each 
leg, i.e. Xj a — > A a , so that the constraint is satisfied only 
on the average on each leg of the ladder. 
Introducing Fourier transformed operators 



fk.a.s y/Tj ^ ' 



b k ,a = -7=YA a e tt 



the mean field Hamiltonian reads 



(14) 



5 



H MF = ( H k + H k) + L \ E {\ J ( A ' + Xl) + U Xa B a ) + ~ J 3 (A» + *§) + 4i' X3 i? 3 



The quadratic terms H\ and Hi are given by 



(15) 



H,'. 



u k,i 
bl, 2 



-4txi cos k — Ai+ (ii -2t'x3 

- 2t'x3 -4tx2 cos k - A 2 + /i 2 




/fe,2,t 
\/-fe,2.J./ 



where and A& are the following 2x2 matrices 



Cfc A fc 
A fe -4 



/fc,2,T 

■'-M.J. 

V/W 



-Ai - ( 2tB x + | J xi ) cos fc -f fl 3 - jJsXs 



t'B. 



A, 



2tB 2 



JX2 COS fc 



A fc 



-§JAicosfc -|JA 3 



JA 3 



t J At COS k 



The mean fields are determined by self-consistently solv- 
ing the single-particle problem of Hmf and calculating 
the corresponding expectation values according to Eqs. 

lima . 

Diagonalization of the bosonic part of the Hamilto- 
nian yields two holon bands. In the ground state of the 
system, the holons are assumed to Bose condense into 
their lowest energy stated Denoting the amplitudes of 
the lowest holon state on the two legs by A a , the bosonic 
bond mean fields become 

B a = 28 A 2 a , a =1,2, (16) 

B 3 = 2SA 1 A 2 . 

The spinon part of the Hamiltonian can be diagonalized 
by a Bogoliubov transformation, from which the self- 
consistent equations for the mean fields A and x> an d 
the Lagrange multipliers A are determined numerically. 
Furthermore we define the BCS order parameters A', as 2 ^ 

A ' 3 ,a: 3 ',a> : = (Cj^fCj'.aM.) (17) 
~ (^],a^j',a')(fj,o-,lfj',a'l) = Bj^ji >a ' Aj,a;j' ,a' 

in terms of the holon and spinon mean fields. For the 
BCS order parameters we use the same bond labeling 
scheme as in Eq. JT^J. 

The redistribution of charge carriers due to the chemi- 
cal potential difference is implemented via the holon de- 
grees of freedom as can be seen in the mean field Hamil- 
tonian. In this way there is no effect at half-filling. Fur- 
thermore, the constraint and the renormalization of the 
coupling constants induces a non-trivial mutual feedback 
for the charge and spin degrees of freedom. 



B. Results 

Within the above mean field description we are able 
to analyze the behavior of the two-leg ladder for differ- 
ent values of the doping concentration, 5, and chemi- 
cal potential difference, Afj,. In particular, we are inter- 
ested in the BCS order parameters as the indication for 
Cooper paring. In the following the parameters of the 
t-J model are fixed to isotropic coupling with t' — t and 
J' = J = 0.5t. Furthermore, we restrict ourselves to the 
low doping region 8 < 0.25, where the above spinon-holon 
decomposition is expected to be qualitatively reliable^. 

For A// = our calculations agree well with the overall 
behavior obtained from similar mean field calculations 
using a Gutzwiller-type renormalization method; 19 While 
at half-filling (8 = 0) the BCS order parameters vanish, 
they increase monotonically with hole doping away from 
half filling. Their values on the legs coincide (A^ = A 2 ), 
whereas a phase shift of ir exists relative to the rung 
order parameter A 3 , in analogy to the d-wave pairing 
symmetry on the square lattice version of the doped t-J- 
model. 

In order to analyze the influence of a chemical po- 
tential difference between the two legs on the Cooper 
pairing, we follow the behavior of the BSC mean fields 
A[ 2 3 upon increasing A[i > for two fixed hole con- 
centrations, 8 — 0.1 and 6 = 0.2, shown in Fig. [31 (a, b). 
In both cases the chemical potential difference leads to 
the reduction and eventual destruction of the BCS mean 
fields. However, there is a qualitative difference between 
the two doping levels. For 8 = 0.1 a crossover from a 
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FIG. 3: BCS mean fields A'(a,b) and hole densities B of the 
t-J model on a two-leg ladder as functions of the chemical 
potential difference A/z, at constant 8 = 0.1 (a,c), and 0.2 
(b,d), for t — t, and J' = J = 0.5t. Values for the lower 
(upper) leg are plotted with solid (dotted-dashed) lines, and 
the BCS mean field on the rungs using dashed lines. 



strong to a weak superconducting regime occurs, while 
no such regime change takes place for S = 0.2. The 
crossover at A/j, w 0.6t in Fig. [3] (a) for 6 — 0.1 co- 
incides with the almost complete hole-depletion of the 
upper leg in Fig. [3] (c). This behavior can be under- 
stood by the following properties of doped two-leg lad- 
ders. These systems constitute typical examples where 
superconductivity originates from a doped RVB phase, 
which is characterized within this mean field approach 
by a finite gap in the spinon spectrum^. Furthermore, 
this spinon gap decreases upon doping holes into a half- 
filled two-leg ladder—. Within the mean field approxi- 
mation we can obtain the spinon gap at finite values of 
A/i for the doping levels considered above. In Fig. 0] the 
development of the spinon gap upon increasing the chem- 
ical potential difference, A/i, is shown for both 6 = 0.1, 
and 6 = 0.2. In cither case does the inbalance in the 
distribution of holes between the two legs lead to an ad- 
ditional reduction of the spinon gap. This behavior is 
expected as the RVB phase in the ladder is dominated 
by the formation of rung singlet pairs. Concentrating the 
holes onto a single leg destroys statistically more rung 
singlets than distributing them equally among both legs. 
However, at 6 — 0.1, a large spinon gap is found even 
for e.g. A/i/i = 2, where the upper leg is almost com- 
pletely depeted, as seen in Fig.|3](c). Although the holes 
are already strongly concentrated onto the lower leg, the 
spinon gap is not destroyed until A/i becomes as large as 
A/t c w 2.6t. Along with the RVB state, d-wave supercon- 
ductivity thus prevails up to this critical value of A/i c , 
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FIG. 4: Spinon gap As of the t-J model on a two-leg ladder as 
a function of the chemical potential difference A/i at constant 
hole doping 5 = 0.1 (solid line) and S = 0.2 (dashed line) for 
t' — t, J' = J — 0.5t. The inset shows the gapless spinon 
bands in the normal state at S — 0.2, Afi/t = 2.0. 
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FIG. 5: Low doping phase diagram of the inhomogeneously 
doped two-leg ladder in the mean field approximation of the 
t-J model at t' = t, J' = J = 0.5t. The BCS mean fields 
vanish beyond the solid line. In the low doping region the 
crossover regime connecting the d-wave SC and a regime of 
reduced BCS mean fields is indicated as a shaded area. 



as seen in Fig. (a). For the larger doping of 8 = 0.2 
however, concentrating the holes onto the lower leg sup- 
presses the spinon gap completely, therebye destroying 
the RVB state. This follows from a comparison of the be- 
havior of the spinon gap in Fig. 0] with the corresponding 
behavior of the charge distribution shown in Fig. 01(d). 
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Along with the spinon gap the superconducting state dis- 
appears already at A/i c f» 1.5t. However, the chemical 
potential difference which is necessary to pull the holes 
onto the lower leg increases upon increasing the hole dop- 
ing level. 

The resulting phase diagram in Fig. OD displays a pe- 
culiar structure. For low-doping concentrations 5 < 0.15 
we observe two regimes, strong and weak superconduc- 
tivity, separated by a broad crossover. The crossover 
region is characterized by the depletion of holes from the 
upper leg. For larger dopings 5 > 0.15 only the regime 
of strong superconductivity remains, and the RVB state 
is destroyed once the holes are sufficiently unequal dis- 
tributed among the two legs. The mean field solution 
suggests that there exists a critical doping S c f=s 0.08 be- 
low which the superconducting state along with the RVB 
spin liquid state remains stable for all A/i. 

Within the mean field approximation the non- 
superconducting phase appears to consist of two indepen- 
dent subsystems. This can be referred from the inset of 
Fig. 01 which displays the gapless spinon bands at 5 = 0.2 
and Afi/t = 2.0, well inside the normal phase. The spec- 
trum consists of the spinon bands of a spin-1/2 antiferro- 
magnetic Heisenberg chain, with nodes at k = ±7r/2, and 
two additional bands, with nodes at kp = ±7r/2(l — 25). 
These nodes correspond to those of a single chain Lut- 
tinger liquid. The system in the normal state therefore 
appear to be separated in a t-J chain with hole doping 25, 
having the properties of a Luttinger liquid, and a spin-1/2 
antifcrromagnetic Heisenberg chain. Furthermore, gap- 
less charge excitations only exist for the Luttinger liquid. 
This complete separation is likely an artefact of the mean 
field approximation, as in the strong rung coupling limit 
J' J the system is obviously a single chain Luttinger 
liquid. This conclusion results also from the analysis of 
the weak-coupling Hubbard model in the next section. 

The mean field description of the t-J ladder is in qual- 
itative agreement with the numerical result of the pre- 
vious chapter. In the following section we will analyze 
the Hubbard model by means of a renormalization group 
treatment in the weak coupling regime, which reflects 
the same basic properties of the inhomogeneously doped 
ladder. 



V. RENORMALIZATION GROUP 

In order to complement the analysis of the t-J model, 
in this section we consider the weakly interacting Hub- 
bard model on a two-leg ladder. A renormalization group 
(RG) treatment supplemented by Abelian bosonization 
allows a detailed analysis of the phase diagram in the 
weakly interacting limit, and a characterization of the 
various phases in terms of the low-energy modes. We fol- 
low an approach established on standard two- and N-leg 
ladder systems (i.e. for A/i = 0)i n i 2 7i 2 ?i 2 9 

In the current case of an inhomogeneously doped two- 
leg ladder we consider the weak repulsive limit, < U <C 



t, t' of the Hubbard model, 

H = - t ^2( C ia.,s C 3+ha,s +h.C.) 

/'^V;.:..^ • h.c.) (18) 

j> 

+^ ^ c j,a,1 c j,aA c j,g,l c j,a,l 

with the same notations as is Sec. I. The quadratic part 
of the Hamiltonian (JTHJ, i.e. Eq. (|T%|) for U = 0, can be 
decoupled via a canonical transformation 

t / lTA/xAP t j l±WD t 

where D = ^4t' 2 + A^ 2 . These rung operators inter- 
polate smoothly between the bonding and anti-bonding 
combinations at Afj, = 0, and the original fermions for 
Afi/t' — > oo, where d)- i s — > ct i s . In momentum space 

two bands corresponding to d\ s (k), i = 1,2 result with 
dispersions 

e 1/2 (k) = -2tcos{k)±^D-fL, (20) 

and a bandwidth At. 

Consider now the effect of the Hubbard interaction 
term in Eq. (fT%|l . When both bands are completely sep- 
arated in energy, only the lower band is filled, and at 
half filling (5 = 0) a band insulator is obtained which 
upon doping (5 > 0) becomes an ordinary spin-1/2 Lut- 
tinger liquid (LL). For D < 4tcos 2 (Tr5), both bands are 
partially filled, and inter-band interaction effects must be 
examined. While this proves difficult in general, progress 
can be made upon considering the weakly interacting 
limit. Since the interest is thus on the low-energy physics, 
the dispersions l|20|l can be linearized around the Fermi 
points kF.i, i = 1,2, determined by £i(kF,i)=£2{kF,2) 
and fcjri + kp^ = 7r(l — 5). Furthermore left- and right 
movers, d^ R ^ Lis , are defined with respect to the Fermi 
level in each band, i = 1,2. For generic (i.e. incommen- 
surate) Fermi momenta the interactions consist of intra- 
and inter-band forward- and Cooper- scattering. These 
can be organized in terms of the U(l) and SU(2) current 
operators 

s 

Jpij = 2 ^""^ < ^p,i,8°'a,a'^P>3i» / i (^1) 
s,s' 

where p = R,L, and the band indices i,j = 1,2. The 
following non-chiral current-current interactions are al- 
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lowed by symmetry 
Hi 



E 



Jdx {c^JbmJlu - cfiJuu ■ Jlu) 

+ E j' Lr Kj J Rij J Lij - <V'/>V. ' Jitj) 



(22) 



Ma 

In this representation / (c) denotes couplings related to 
forward- (Cooper-) scattering, and the symmetry of the 
inter-band scattering terms under the band exchange is 
explicitly taken into account. Using current algebra and 
operator product expansions, a one-loop RG flow for the 
various couplings can be derived^i which in our notation 
reads 



dl 

dl 
dc p 12 
dl 



dcf 2 
~df 



dl 
dl 



1 

2^ 
1 

2m 



( C P)2 + Ji( c -.)2 
v ix> 16 



2 

E- 

^ 2v 



v\ + v 2 

E- 

4- 2v t 



ir' 7 ) 2 + 2c" 



er\2 



2u, 



■(cS) 



^ r P _|_ —r- 17 ^ 
C \i C i2 ' ^Q C li c i2 



16 



-li c i2 



1 



' C li C i2 + 2 CliCi2 



(23) 



vi ■ 
1 



(-2 



2 T c 12/l2 2 C< ^ 2 ^ 2 



1 



i'2 



(c 12 ) 2 + ^(c 12 ) 2 

2 C 12 C 12 — o( c l 2 ) 



(/r 2 ) 2 



where i = 1, 2, 1 = 2, 2 = 1, and = 2isin(fcp.i) are 
the Fermi velocities for the bands. The successive elimi- 
nation of high frequency modes is obtained from l|23l) by 
integration along the logarithmic length scale related 
to an energy scale E ~ te~^ 1 . The flow equations can 
be integrated once the bare values of the couplings are 
known. For the Hubbard interaction of Eq. I|18|) they are 
obtained as 



-li 



L 22 



- "° = 4c? x = 4c p 22 = U 



1 + 



<2 = ft 



4/" 



12 



u 



Ap 

IT 

Afj, 
If 



(24) 



Increasing Ap away from zero can be seen to reduce 
the bare inter-band scattering with respect to intra-band 
scattering. 

Depending on the parameters, integration of the flow 
equations leads to different asymptotic behavior, with 



either a flow to a finite-valued fix point, or to instabili- 
ties characterized by universal ratios of the renormalized 
couplings beyond a scale I*, where the most diverging 
coupling becomes of the order the bandwidth. While the 
consistency of the one-loop renormalization group equa- 
tions is restricted to I < I*, the asymptotic ratios can be 
utilized to derive a description of the low-energy physics 
of the system within Abelian bosonization^S Introduc- 
ing canonically conjugated bosonic fields and W vl 
for the charge and spin degrees of freedom [y = p, a) 
on each band i = 1,2, the fermionic operators can be 
represented as 



1 R/L,i,s 



V 



is_ „iy/^/2[±($ l ,i+s<S>„ i )-(9 pi +sO„i)] 



V2i 



(25) 



where d V i is the dual field of $,/,, so that d x 9 u i = n„j. 
The r]i s are Klein factors, ensuring anticommutation rela- 
tions, and a is a short-distance cutoff^i Using the above 
representation, the interacting fermionic Hamiltonian 
transforms into a bosonic Hamiltonian, Hb = H q + Hi, 
containing quadratic terms 



H q = £- dx 



Pv\Pv 



dx 



fl2 



-{d x <$>vid x $v2-Ti vl Il v2 ) 



PAPA 

and sine-Gordon-likc interaction terms 



H 



(26) 



Hr = 



Jdx{ cf x cos (V2l3 a $ ai ) + c a 22 cos {V2f3 a <S> a2 ) 

-4c? 2 cos (2/? p 6> p _)[cos - cos {fi a 6 a J)\ 

-cf a cos (2/3 p p _)[2cos + cos 
-c a 12 cos (2/3 p 6» p _) cos (/30 ff _) 

+2/f 2 cos (AA-) cos }, 



(27) 



where (3 p — y^r, fi a = — V4"7r, and the fields $j,± = 
($„i ± *i/a)/V2 and n„± = (n„i ± tt v2 )/V2 have been 
introduced. Upon minimizing the energy in a semiclassi- 
cal approximation, any coupling that diverges under the 
RG flow opens up a gap for a field that is pinned by the 
corresponding terms in (|27|l . 

Performing the above procedure, four different phases 
are obtained for the Hamiltonian of Eq. i|18|) . shown 
in the (<5, A/i)-plane for isotropic hopping, if — t, in 
Fig. 5. The various phases are labeled according to the 
number of gapless charge (n) and spin (m) modes by 
CnSm. The different asymptotic regimes of the RG flow 
(|23|l are related to the phases shown in Fig. 5 as follows: 

(C1S1) The single band LL. This phase with an empty 
upper band is labeled C1S1, reflecting the number of gap- 
less modes. For incommensurate filling the dominant 
correlations are charge density waves (CDW) and spin 
density waves (SDW)i2i 

(C2S2) The trivial fixed point. In this regime the cou- 
plings stay of order U under the RG flow, or renormalize 
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FIG. 6: Phase diagram of the inhomogeneously doped Hub- 
bard model on two-leg ladder in the weakly interacting limit, 
and for t' = t. The phases CnSm are labeled according to the 
number of gapless charge (n) and spin (m) modes. Solid lines 
are phase boundaries. The dashed line indicates the crossover 
from a region with dominant SC correlations in the lower and 
CDW and SDW correlations in the upper band (C2S2.I), to 
a region dominated by CDW and SDW correlations in both 
bands, for larger values of A/i (C2S2.II). 



to zero. Therefore no gap opens in this phase, labeled 
C2S2. Furthermore, the two partially filled bands are 
decoupled in this regime. Standard LL theory«2i, used to 
determine the dominant correlations within each band, 
indicates a crossover between two regions labeled C2S2.I, 
and C2S2.II respectively, c.f. Fig. 5. The dominant 
correlations in the C2S2.II regime are CDW and SDW 
within both bands, whereas in the C2S2.I regime the 
lower band is dominated by superconducting (SC) fluc- 
tuations. 

(C2S1) Single-band superconductivity. Here, all cou- 
plings stay of the order of U or renormalize to zero, except 
for ~ ~V2- This results in a pinning of the spin mode 
of the lower band, and the number of gapless modes re- 
duces to C2S1. SC correlations dominate the lower band, 
and CDW and SDW correlations the upper band. Fur- 
thermore, inter-band phase coherence is not established 
within this regime. 

(C1S0) Inter-band superconductivity. In this regime 
the diverging couplings flow towards the asymptotic ra- 
tios 



gapless and exhibits dominant superconducting pairing 
correlations, with a sign difference between the bands. 
This is usually referred to as the d-wave -like supercon- 
ducting phase of the two-leg ladder / 1 ' 27 

The phase diagram in Fig. 5 confirms the results 
obtained along the line A/i = Om^^ But it also in- 
dicates a limited stability of the various phases found 
at A/i = under inhomogeneously doping of the lad- 
der. Upon increasing A/i, superconductivity is gradu- 
ally suppressed, with intermediate phases showing resid- 
ual superconducting fluctuations. Consider for example 
the low doping region where inter-band d-wave super- 
conductivity occurs for vanishing A/i. While for small 
A/i > d-wave superconductivity sustains, inter-band 
phase coherence is lost when A/i reaches a value of ap- 
proximately 1.51/ . For larger values of A/i intra-band 
superconductivity persists within the lower band (which 
predominately projects onto the lower leg of the ladder). 
In the upper band spin fluctuations have become gap- 
less and SDW and CDW correlations dominate. Further 
increase of A/i results in the suppression of all supercon- 
ducting correlations, giving rise to two-band LL behav- 
ior. In the weakly interacting limit the upper band is 
depleted for Aii/t > 3.5, and a single-band LL, resid- 
ing mainly on the lower leg of the ladder, dominates the 
large-A/i regime. This progressive reduction of super- 
conducting pairing correlations is also observed in the 
finite-temperature phase diagram. In Fig. 6 we show re- 
sults for t' = t, and Afi/t = 0.3 in the (5, T)-plane. From 
the renormalization of the energy scale, E ~ te~ nl , the 
logarithmic length scale of Eq. i|23|) can be related to a 
temperature scale T = E = Toe~ 7ri . While the phases 
of the system for T — > are found in accordance with 
Fig. 5, the finite temperature phase diagram reveals a 
successive enhancement of superconducting pairing cor- 
relations with decreasing temperature. Consider again 
the behavior close to half filling. At high temperatures 
the system is dominated by LL behavior in both bands. 
Upon decreasing the temperature gapless superconduct- 
ing correlations develop within the lower band. At even 
lower temperatures, a finite spin gap opens for the lower 
band, then finally phase coherent inter-band d-wave su- 
perconductivity emerges, along with the opening of the 
second spin gap. Thus in an intermediate temperature 
regime, well above the onset of d-wave superconductiv- 
ity, a single spin gap persists in the lower band, which is 
related to the bonding band at small values of A/i. This 
partial spin gap formation might be interpreted as a phe- 
nomenon similar to the pseudogap phase in the HTCS 
materials. 



4^ = 8/^ = ^2, clJvx = c^/ua. (28) 

From the low-energy effective bosonic Hamiltonian two 
finite spin gaps are obtained, and a pinning of the charge 
mode Q p - = 0, resulting in phase coherence between the 
two bands. The number of gapless modes is reduced to 
C1S0. The remaining total charge mode, (& p+ ,8 p+ ), is 



VI. DISCUSSION AND CONCLUSION 

For many years two-leg ladder systems have been 
prominent model systems for discussing superconductiv- 
ity in doped spin liquids. In this paper we investigated 
the stability of the superconducting phase for inhomoge- 
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FIG. 7: Finite Temperature phase diagram of the inhomogc- 
neously doped Hubbard model on the two-leg ladder in the 
weakly interacting limit for t' = t, and A/i/t = 0.3. Solid 
lines are phase boundaries, whereas the dashed line indicates 
a crossover inside the gapless regime. 



neous doping by various approaches which yield an over- 
all similar picture. 

As anticipated from a strong coupling point of view a 
chemical potential difference between the legs of the lad- 
der acts pair breaking, as could be clearly demonstrated 
in numerical exact diagonalization of finite systems. 

The mean field analysis based on the spinon-holon de- 
composition suggests that the inbalanced carrier distri- 
bution indeed leads to the suppression of the supercon- 



ducting state on the doped ladder. Nevertheless, a more 
differeniated picture emerges. In the low-doping region 
the RVB state remains stable even for large differences 
in the chemical potential and supports a weakly super- 
conducting phase. This RVB phase and the weak super- 
conducting state do not exist for higher doping concen- 
trations above 6 w 0.15. 

A modified picture is observed in the renormaliza- 
tion group treatment of the weakly interacting Hubbard 
model on the two-leg ladder. Also here inhomogenous 
doping leads to a suppression of the superconducting 
phase, a Luther-Emery-liquid characterized by one gap- 
less charge mode (C1S0). Moreover, an intermediate 
phase appears which corresponds to a single channel be- 
ing superconducting while a coexisting channel forms a 
Luttinger liquid (C2S1). In both the t-J and the Hub- 
bard model a phase of complete destruction of supercon- 
ducting fluctuations appears for large enough differences 
in the chemical potential. Within the renormalization 
group approach this normal phase is characterized as a 
single Luttinger liquid state (C1S1). While this identifies 
the true low-energy properties of this regime, the change 
in the spinon spectrum discussed in Sec. Ill rather re- 
flects a short-coming of the mean-field solution. 

In conclusion we emphasize that inhomogeneous dop- 
ing of the two-leg ladder is harmful for the formation 
of the superconducting state. Furthermore, it can be an 
interesting tool to access new phases for this type of elec- 
tronic systems. 
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Haas, Masashige Matsumoto, and Igor Milat. This work 
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fund. 
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